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Abstract 

A primordial magnetic field may be generated during an inflationary pe- 
riod if conformal invariance is broken. We reexamine and generalize previous 
results about the magnetic field produced by couplings of the form R n F tlv F^ v . 
We show that the amplitude of the magnetic field depends strongly on n. For 
adequate values of n the field produced can serve as seed for galactic magnetic 
fields. We also compute the effective interaction between the electromagnetic 
field and the geometry in the context of scalar QED (with and without classi- 
cal conformal invariance). In both cases, the amplitude of the magnetic field 
is too small to be of astrophysical interest. 
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I. INTRODUCTION 



Magnetic fields play an important role in a variety of astrophysical situations. There is 
enough evidence for the existence of intragalactic magnetic fields , with an amplitude of 
10~ 6 G and uniform on a scale of lOkpc. It is not completely clear how these magnetic fields 
were generated. A plausible explanation is that some kind of dynamo effect || could have 
amplified a pre-existent magnetic field. But then the question is about the mechanism that 
produced this 'seed' field. 

An attractive suggestion is that it has a primordial origin and could have been produced 
in the early universe during an inflationary period @-|6|]. Denoting by r the energy density 
of the magnetic field relative to the energy density of the cosmic microwave background 
radiation, r = a pregalactic magnetic field caracterized by r ~ 10 -34 is needed in order 
to explain the present value of r ~ 1. 

As pointed out by Turner and Widrow H, it is not possible to produce the required 
seed field with the usual Maxwell Lagrangian. The reason is conformal invariance. Indeed, 
in a conformally invariant theory B decreases as \, where a is the scale factor of the 
Robert son- Walker metric. During (exponential) inflation, the total energy density in the 
universe is constant, and the magnetic field energy density is strongly suppressed, giving 
r(A) = 10- 104 (1M2£) 4 . 

Conformal invariance can be broken in different ways. From a phenomenological point 
of view, one can consider Q couplings of the form RA^A^. This type of interaction-terms 
give rise to the required seed, but, not being gauge- invariant, are theoretically unappealing. 
In string-inspired models, conformal invariance is broken by the coupling between the elec- 
tromagnetic field and the dilaton. This coupling may produce the seed field ||. Finally, 
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one can consider gauge-invariant couplings of the form ^RF llv F tJ,v . However, the seed field 
produced is extremely small [|J. 

In this paper we will reexamine the generation of primordial magnetic field due to the 
above mentioned gauge-invariant couplings. These terms appear due to quantum effects 
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when taking into account one loop corrections for QED in curved spaces |7[ (throughout 
the paper we will work in the context of scalar QED). In an expansion in powers of ^ 
(Schwinger DeWitt Expansion, SDWE ||) [], one expects the effective action to contain 
couplings the form -^F^F^ . But as already mentioned in Ref. [[|, during the inflationary 
period one tipically has TZ ^> m 2 . Therefore, there is no reason to keep the lowest order 
contribution n = 1, and it is of interest to investigate more general couplings. 

In the next Section we will compute the primordial magnetic field produced by couplings 
of the form ^F^F^. For n = 1, our results are considerably smaller than those of Ref. 
Q. The discrepancy is due to an overestimation of the value of r at first horizon crossing. 
However, we will show that the amplitude of the magnetic field depends dramatically on n. 
In particular, for adequate values of n, it is possible to generate a sufficiently large seed field 
to explain the present values of the galactic field through a dynamo mechanism. 

In this situation, the obvious question is about the effective Lagrangian in the opposite 
regime TZ 3> m 2 . We address this issue in Section 3. Using an improved version of the SDWE 
ITJ , we will show that, in the leading-logarithm approximation, the effective coupling is of 
the form F^F^ In We will also compute the amplitude of the magnetic field produced by 
this coupling. Unfortunately, since the logarithm is a slowly varying function, this amplitude 
will be extremely small. 

In Section 4 we discuss the case of QED with conformally- invariant quantum fields ||. 
Again, we will find that the value of r is too small to be the seed of the galactic dynamo. 

Throughout the paper we will use units in which h = c = 1. 

II. PRIMORDIAL FIELDS AND THE SCHWINGER DEWITT EXPANSION 

The Lagrangian for scalar QED on curved backgrounds is given by 



We denote by TZ any component of the Riemann tensor 
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£ = —F^F^ - D^* + m 2 00* + £ R #* (1) 

where _D M = <9 M — ieA^ is the covariant derivative for the scalar field. The theory is confor- 
mally invariant only for massless and conformally coupled (£ = §) fields. 

One can compute an effective Lagrangian for the electromagnetic field by integrating out 
the quantum scalar field. Using dimensional regularization, the effective Lagrangian can be 
expanded as follows U 

1 1 1 /m\ d " 4 00 d 

£e ff = -± V" + ^ E°i(^) ™ 4 " 2j r (j - §) (2) 

4 2 (47r) 2 \ n J j=Q I 

where the a,j(x) are the SDW coefficients, and d is the spacetime-dimension. The first three 
terms are divergent in the limit d — *■ 4, and the poles must be absorbed into the bare 
constants of the classical Lagrangian and into a redefinition of the electromagnetic field. 
The first SDW coefficients for a charged scalar field are given by fTO 



a = 1 

ai = -(£-£) 

d2 = igo [R^vaP R^ va ^ — R^R^^j + 

rto — _ e 2 f fp^p _ e 2 fp^ p ;p 
— £2 (nF ] F pv 4- <= 3 F F up F^ 

30 V ~ 30 M 1 ' 1 P 

_|_ e 2 D PPO" _ e 2 D FPP 

^ 60 P^P°" 1 90 P y 1 P 

+^RF, U F^ - ^RF^F^ + ■■■ (3) 

where we omitted the purely gravitational terms in 03. 

After absorbing the poles, the effective Lagrangian contains the usual Maxwell term plus 
the finite corrections 



m 2 2 m 2 m 2 a fe 

a Q m ln— + a im In — + a 2 In — + ^ — — — 2 T {k - 2) (4) 

P P A* fc>3 I'"' 
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From the above equations we see that the first non trivial interactions between the gravita- 
tional and electromagnetic fields are contained in a 3 and read 

RF^ , R^F^F" , R^F\ (5) 

In the same fashion, the SDW coefficient a n+2 contains interactions that can be (schemati- 
cally) written as ^^F 2 . 

In what follows, we will analyze the magnetic field produced by a typical term in the 
SDWE. For simplicity we consider the following effective Lagrangian 



4// = ~^(i + l^f] (6) 



The classical equations of motion are given by 



+ b (-§)") = (7) 

fit 

In the particular case of a spatially flat Robertson Walker metric 

ds 2 = a 2 (rj)(-dr} 2 + dx 2 + dy 2 + dz 2 ) (8) 

they read 

&*(F IW (1 + b(— 2 D = (9) 

Obviously, if R is constant (exponential inflation), these equations do not differ from Maxwell 
equations. Therefore, non trivial effects appear only for extended inflation. 

We will work in the radiation gauge A = J2i=i diAi = The Fourier expansion of the 
field is 

r d 3 k 2 

MV,*) = J (2^X>f ® [a W (k)A k (r))e^+ h.c. ] (10) 

where u = k, a^(k) and a^ x \k) are the usual annihilation and creation operators. The 
vectors e ( x \k) satisfy e ^(k). k — 0. The classical equation for the Fourier modes reads 

(A k + k 2 A k ){\ + b (— 2 T) + -^-R^A k R = (11) 
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where the dots denote derivatives with respect to the conformal time rj. 

During the inflationary period, R 3> m 2 and the quantum correction dominates over the 
Maxwell Lagrangian. Indeed, from Einstein equations we have R ~ 0(/?t ot /m^) where p to t 
is the total energy density. Consequently, the term b(R/m 2 ) n F flu F t11 ' dominates over the 
usual F ilv F llv as long as p tot ^> m^m 2 . This gives p tot ^> (10 8 GeV) 4 for m = m e , 
the electron mass. In this situation, Eq. |ll] reads 



R 

A k + k 2 A k + n-A k = (12) 
ri 

Before proceeding, we would like to stress an important point related to the normaliza- 
tion of the Fourier modes A k . As our effective Lagrangian is now —\{R/ 'm 2 ) n F^ U F^ U , the 
canonical commutation relations between the coordinates and the conjugate momenta 
n * = A^bC^Y read 

R \ n ~ 



(R \ 
m 2 ) 



iUi-^i) S 3 (x-x') (13) 



The usual commutation relations between the creation and annihilation operators are com- 
patible with the above commutator only if 

A k A* k -A k A* k = 4^r (14) 

This will be important in what follows. 
Consider the quantity 

<^8^< B2 > < 15 > 

where (• • ■) denotes the vacuum expectation value. Although this is not the magnetic field 
energy density associated with the effective Lagrangian (|6|), it becomes the energy density 
when C^p) < 1 As a consequence, we can analyze the time evolution of in order to 
obtain the value of r at a time when the Maxwell Lagrangian dominates. For simplicity, we 
will refer to ps as the magnetic energy density. 
In terms of the Fourier modes we have 
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Pb = ^~ 4 f-^-2k 2 \A k ( V )\ 2 (16) 



\na A J (2tt) 3 2o; 

The energy density in the fcth mode of the magnetic field, defined as ps(fc) = kcLps/dk, is 
given by 

PB(k) = — 1-j A: 4 lAkW I 2 (17) 
16 7T d a 4 

During extended inflation we can think the Universe as filled with a perfect fluid having 
an equation of state p = jp, with —1 < 7 < —1/3. The scale factor evolves as a(t) = 
where a = 2/3(1 + 7) > 1 and T is chosen so that atoday — 1- This implies that, in conformal 
time, R/R = 2/(a — 1)77. For this evolution, the properly normalized solution of Eq.[l2| is 
given by 

n_ n 

AM - (-*,)- f r ^n ~ (^) /#><-*,) («) 

where v = | — and H^\x) is a Hankel function. In order to select a unique solution 
we have also imposed that A k (rj) ~ e _lfcr? as 77 — > —00. This is the natural choice for the 
vacuum state in the in region. 

First horizon crossing takes place when — kn = 2ixa/(a — 1). Inserting this value of n in 



Eqns.|17| and |TJ| one finds the value of pb at that time. The value of ptot(v) can t> e obtained 



from Einstein equations at first horizon crossing. It is given by 



2 . 2 ■>„ , . _2 



m z pl /47r\«-i /3a\"-i / A 



a— 1 



«-w=6?lfJ ItJ sf T (19) 



The value of r at first horizon crossing, denoted by r , is therefore given by r = 



After first horizon crossing the quantum fluctuations are assumed to become classical per- 



turbations. For fci; > 1 the term k 2 A k in Eq.[12] can be neglected. The time evolution of the 
Fourier modes is approximately given by 



Ti S ' 

A k + A k = (20) 

77 



where s = —6(1 + 7)/(l + 37). 



From here, the analysis follows closely that of Ref. Eq.^ admits a constant solution, 
which gives an uninteresting p B ~ a~ 4 . It also has the solution 

A k ~ J^lt^*™ ~ ( a )^+i)+(en+3 )7] (21) 
The solution grows with time for 7 > — |^r|- In this situation, and as long as R 3> m 2 , 

2n + 27(n + l) 

r evolves like p tot 1+7 . When ~ m 2 , or when conductivity effects during reheating 

(1-37) 

0| become dominant, ps ~ a -4 and r oc Ptot + ^ ■ Assuming a very rapid reheating, the 
Universe enters the radiation dominated period, ptot equals p 7 , and r becomes a constant. In 
order to calculate this final value of r we have to distinguish between two possibilities. The 
first one corresponds to the case in which the growing solution disappears before the plasma 
effects become dominant. This means that M 4 < b~ '*» m 2 ; m 2 , where M 4 is the total energy 
density at the end of the inflationary period, and b~™ m^m 2 is the total energy density at 
the instant in which the R n F 2 term ceases to dominate over the usual Maxwell term. In 
this case, and using the relations betwen r and p tot stated above, we get that the value of r 
at the beginning of the radiation dominated period is given by 

2n + 2(n + l)7 1-37 

\b n m vl J \b ™ m vl m z I 

The other possibility is that the conductivity becomes important before the growing solution 

2n+27(n + l) 

disappears, i.e., that M 4 > b~~m 2 l m 2 . Then, r oc p tot 1+7 during all the period that 
goes from first horizon crossing to the end of inflation, and therefore the final value of r 
reads now 



' ptot{k)' 



2n + 2(n + l)7 



r = r„^j (23) 

For b = 0(1) and m = m e , we found that there is a value of 7 between —0.7 and —0.6 that 
separates the two cases. The results for different values of 7, A = 1 Mpc and m = m e are 
summarized in Table I. 
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Table I 
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r 


1 Mpc 




n - 


= 1 


n = 2 


n = 3 


n = 4 


n = 5 


-0.5 


10" 


-84 


1Q -67 


IO" 51 


1Q -35 


io- 19 


-0.6 


10' 


-92 


1Q -63 


10~ 37 


io- 10 


10 17 


-0.7 


10" 


116 


1Q -96 


io- 77 


1Q -56 


1Q -36 


-0.8 


(io- 


-141\ 


1Q -146 


10 -147 


1Q -148 


10 -149 


-0.9 


(io- 


-184 \ 


(io- 222 ) 


10 -253 


10 -269 


10 -285 



Table I: r\iMpc as a function of 7 and n. The numbers between brackets correspond to 
situations where there is no growing solution. The value of r for n = 5 and 7 = —0.6 cannot 
be trusted because it does not satisfy the condition r < 1. 



These results depends strongly on n. For 7 > — .6, we are in the case in which the growing 
solution disappears before the conductivity becomes dominant. We can see from Table I 
and from Eq.^ that the value of r increases exponentially with n (for n high enough, the 
resultant field can serve to seed the galactic dynamo). For 7 < —.7, plasma effects dominates 
before the growing solution disappears, and therefore the period of time during which the 
amplification mechanism is operative, is less than the corresponding one in the opposite 
case. Then, the value of r as given by Eq.[23|, is a compromise between the amplification and 
the time during which it takes place. In Table I we can see this behavior: for 7 = — .7 we 
still have an exponential growth of r, while for 7 < — .8 the period of amplification is not 
enough and r now decreases exponentially with n. 

It may seem strange that we obtain an exponential decrease in certain cases while we 



can see from Eqs.p2| and |23] that the amplification factor grows exponentially with n. The 
cause of this behavior is in the dependence of the initial condition r on n (see Table II). 
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Table II 


7 


r o 1 Mpc 




n = 


= 1 


n - 


= 2 


n = 3 


n = 


= 4 


n 


= 5 


-0.5 


10 


-41 


10 


-58 


1Q -76 


10' 


-94 


10 


-112 


-0.6 


10 


-43 


10 


-71 


10 -ioo 


10" 


-128 


10 


-157 


-0.7 


10 


-44 


10 


-78 


10 -112 


10" 


-146 


10 


-180 


-0.8 


10 


-43 


10 


-78 


10 -113 


10" 


-148 


10 


-183 


-0.9 


10 


-45 


10 


-83 


10 -121 


10" 


-159 


10 


-197 



Table II : 

r o|i Mpc for the electromagnetic field, for different values of 7 and n 

As we see, r decreases exponentially with n, so the behavior of the final value of r with n is a 
compromise between the decrease of the initial condition and the growth of the amplification 
factor, both of which depend on 7. 

At this point, we would like to compare our results for n = 1, with those of Ref. [|J. Our 
values for r are smaller by several orders of magnitude (compare our r{p/ = —.6) ~ 10 -92 
with the correspondent result r ~ 10~ 68 in Ref. (3J). ^ i s easy to find the origin of the 
discrepancy. The calculation in Ref. 0] begins after first horizon crossing and assumes that 
the initial value of pB is H 4 . This is the value of the energy density for a massless minimally 
coupled scalar field. Using this value the resulting r varies between 10 -33 and 10 -12 (for 
—0.9 < 7 < —0.5). Here we computed PB\ist crossing from first principles, and found much 
smaller values (see Table II). The physical origin of the discrepancy is the normalization 
condition Eq.|i~4]. 

Up to here, we have made our calculations using m = m e . Now we turn to study the 
dependence on m of our results. The mass m enters the calculation at two different points. 
First, it determines when the F 2 coupling ceases to dominate over the usual Maxwell 

term, and hence it contributes to determine the amount of amplification that will take place. 
Second, it appears in the initial condition in such a way that r oc m 2n . We can see that 
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raising the value of m improves the initial condition r . On the other hand, it makes shorter 
the amplification period, and the final result is a compromise between these two effects. We 
have analysed the mass dependence taking these facts into account. We have found that 
for M 2 < rripim, the results decrease with increasing mass, while for M 2 > m p im 
the results increase with mass. In the first case the shortening of the amplification period 
preponderates over the increasing of the initial condition, making the result decreasing with 
increasing mass. In the other case, the duration of the amplification period depends only on 
the value of M and no longer on m. This means that raising the value of the mass improves 
the initial condition but does not affect the amplification, giving an increasing result with 
increasing mass. However, in this case the value of m cannot be raised arbitrarily, because 
it must obey the relation M 2 > b~« m p im. 

To summarize, although the value of r is strongly suppressed with respect to that of 
a minimally coupled scalar field, and although the suppression increases with n, there are 
values of 7 such that, after first horizon crossing, a large amplification takes place, and 
produces final values of r large enough to serve as the seed magnetic field. 



III. IMPROVING THE SDWE KILLING THE MAGNETIC FIELD 

From a phenomenological point of view, the couplings discussed in the previous section 
may help to solve the problem of the generation of a primordial field whitout breaking gauge 
invariance. Do they have a theoretical motivation? In order to answer this question, we 
need to compute the effective action for the electromagnetic field in the limit 1Z 3> m 2 . To 
do this, we will use the improved version of the SDWE developed by Parker and Toms in 
Ref. [|TI]. 



A partial ressumation of the SDWE can be achieved by doing the expansion in inverse 
powers of rh 2 = m 2 + (£ — \)R instead of inverse powers of m 2 . In this conjectured 
in Ref. [IT| and proved in Ref. the new SDW coefficients (denoted here by bj) do not 



contain powers of the scalar curvature R. The first bj are given by 
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12 

/,„ — _ e 2 F fp^;p _ e 2 f?A tl/ f ;p 

i£D FP^ pp°- _ <£d FPP/7^ + ... ('241 

T 60 M"^ 90 P^ P ' V / 

where we did not include the purely gravitational terms in 6 3 . Note that &i and bi do not 
contain the terms proportional to R and R 2 that are present in a\ and 02 (see Eq|J). The 
same thing happens for all bj. 

The improved expansion can be obtained from the usual SDWE (Eq.Q) substituting m 2 
by in 2 and the coefficients a,j by the bj. In the limit (£ — ^> m 2 , the coefficient 62 induces 
a coupling of the form 

1 / ^F^U*) (25) 



4\ 967r 2 / \/i 2 



Under the assumption In 3> 1, this term dominates and the effective lagrangian reduces 

p 

to 

where /1 is an arbitrary scale. It is worth noting that this result is also valid beyond one-loop, 
in the leading-logarithm approximation. A similar effective lagrangian has been proposed 
for non-abelian gauge theories in Ref. [|L| . 



There is a simple physical interpretation of the effective Lagrangian given by Eq. |2(| 
If we couple C e ff to an external current, C e ff — > £ e ff + eJ^A^, after a rescaling of we 
obtain the usual Maxwell Lagrangian with a running (curvature-dependent) electric charge 

e2{R > = ,_ ¥ib (27) 

1 96tt 2 111 /I 2 



This is the scale dependence dictated by the renormalization group. This interpretation 

eff 

12 



also makes explicit the fact that C e r t is valid only under the assumption In <C 1 . 



Otherwise, an analysis of strongly coupled QED is needed (strictly speaking, for such large 
values of the curvature one should replace QED by a GUT). 

We will perform the same calculation we did in the previous Section, that is, calculate 
the value of r for a scale of 1 Mpc. Typically, this scale crosses outside the horizon not long 
after curvature falls below the GUT scale (10 28 GeV 2 ). Thus, we cannot use the effective 



lagrangian fl26|) to calculate the initial condition, so we will consider r as given and study 
its evolution after first horizon crossing. 
The equation for Fourier modes is 

A k + k 2 A k %- 7 4^ r = (28) 

We will take p 2 = Rrh, the value of the scalar curvature at the end of reheating, and we 
will consider that e, the electric charge when R = Rrh, is of the same order of magni- 
tude of its value today. Then we can assure that, during most of inflation, In ^> 1 
and g|^2 In <C 1, the conditions needed for the validity of the leading logarithm ap- 
proximation. In this situation, for modes outside the horizon, we can approximate Eq. [28| 
by 

e 2 R ■ 

A k = — A k = (29) 

k 96tt 2 R y ' 

This equation has the same form of Eq. |2(], if we make the identification n <-> — gf^2 = %■ 

The solutions are the same as those of Eq. |20], that is, A k = constant, which gives ps oc a -4 , 

and 

A k ~ ( a )5[( 6 -+ 1 )+( 6 -+ 3 )7l (30) 
Using that x ~ 10~ 4 we can write Eq. |3D] as 

A k ~(ay* [1+ ^ ] (31) 

Since 1 + 37 < 0, this means that, for this solution, pB decreases more rapidly than a -4 . 

So we can consider only the contribution of the constant solution and conclude that for the 

(1-37) 

effective lagrangian given by Eq. ^6], ps decreases as a -4 , and then r oc piot +l) ■ 
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The final value of r will be 

/ fl/f 4 \ 3(1+7) 

(32) 



r = r n 



° \Ptat(k)J 

For 7 < —0.5, the amplification factor that accompanies r is smaller than 10~ 49 . Since the 
initial condition is tipically much less than 1, the value of r is in all cases very far from the 
required value 1CT 34 . 

To conclude, although the higher order terms in the SDWE TZ n F 2 generate large ampli- 
tudes for the magnetic field, once these terms are ressumed they generate only a logarithmic 
interaction InR F 2 . This interaction produces a disco uragingly small seed field. 

IV. TRACE ANOMALY AND PRIMORDIAL FIELDS 

We end the paper with a comment about the conformally invariant case m 2 = 0, £ = 1/6. 
Of course in this situation neither the SDWE nor its improved version are useful to compute 
the effective Lagrangian. However, it is easy to find a closed expression for it. 

It is well known that conformal invariance is broken by quantum effects. The (anomalous) 
trace of the energy momentum tensor is given by the second SDW coefficient || 

r »" = -!^ = ^V + - (33) 

where the dots denote purely gravitational terms proportional to TZ 2 and VV7£. 

In a Robertson Walker space-time, the above trace anomaly can be derived from the 
effective Lagrangian jll] 

C eff = -~F„ v F^(l + K\n^J+.... (34) 

where k = ^2 and a is a reference value for the scale factor (for example the value of a 
after inflation). Here the dots denote non-local terms independent of the scale factor. 



2 There is a sign difference between the logarithmic term of the effective lagrangian in Ref. [14] 
and that of Ref. f§. We agree with the result of Ref. |14| . 
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As before, Eq.|34| can be interpreted as the usual Maxwell Lagrangian with a scale de- 
pendent electric charge 



6 (a) = (l + .ln^) (35) 

V an / 



The effective Lagrangian is valid as long ln^ < 1. 

Under this assumption, one can estimate the magnetic field generated by the conformal 
anomaly. The resulting value for r is again extremely small, of order 10 -104 for A ~ IMpc. 

It has been argued |J that this could be an efficient mechanism if a large number of 
massless (m <C H) fields are present during inflation, so that k — > Nk ~ 0(1). Using the 
approximation 

1 / n \ Nk 

C, tl K --F^ (_) (36) 

it is possible to show || that the previous result for r is multiplied by a huge amplification 
factor A = (f ^ ~ 10 507Vk for F = 10 12 GeK and A = IMpc. 



However, during that period In ^ varies between —60 and and the approximation (|36] 



breaks down for Nk ~ 1. In this situation, even the effective Lagrangian (p^ ) is inadequate, 
since at some point 1 + Nk In — vanishes and one reaches the Landau singularity. Therefore, 
the conclusion is that, if Nk is sufficiently small in order to assure that iVfidn ^ < 1 during 
the period of interest, the amplification factor is very small, A ~ 1 + Nk In ^ . The amplitude 
of the magnetic field generated is too small to be of astrophysical interest. On the other 
hand, if Nk ~ 0(1), the final result for r would depend on the physics of QED at strong 
coupling. 
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